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ABSTRACT 



We present photometry of the Anomalous X-ray pulsar IE 1048.1—5937 in 
the infrared and optical, taken at Magellan and the VLT. The object is detected 
in the I, J and K s bands under excellent conditions. We find that the source has 
varied greatly in its infrared brightness and present these new magnitudes. No 
correlation is found between the infrared flux and spin-down rate, but the infrared 
flux and X-ray flux may be anti-correlated. Assuming nominal reddening values, 
the resultant spectral energy distribution is found to be inconsistent with the 
only other AXP SED available (for 4U 0142+61). We consider the effect of the 
uncertainty in the reddening to the source on its SED. We find that although 
both the X-ray and infrared fluxes have varied greatly for this source, the most 
recent flux ratio is remarkably consistent with what is is found for other AXPs. 
Finally, we discuss the implications of our findings in the context of the magnetar 
model. 

Subject headings: pulsars: individual (IE 1048.1-5937) 



1. Introduction 



Over the last decade or so, evidence has been mounting for the existence of a class 
of neutron stars called the magnetars (Thompson & Duncan, 1996). They have enormous 
external magnetic fields of order ~ 10 14 G (and even larger internal field), the decay of which 
powers luminous high energy radiation. Examples of this type of astrophysical source are the 
Soft Gamma-ray Repeaters (SGRs), which give sporadic bursts of hard X-ray/soft gamma 
rays as well as rare, very luminous (~ 10 44 erg) "giant flares"; and the Anomalous X-ray 
Pulsars (AXPs), so called because their luminosity far exceeds the spin-down value, and no 
binary companion is seen. This implies that the source of power cannot be rotational energy 
and that accretion is excluded. Both the SGRs and AXPs have spin periods of order P = 10s 
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and derivatives P = — 10~ 10 ... — 10 -12 , and are inferred to be young from their energetics 
and spin-down, and from supernova remnant associations in some cases. The two groups 
have been linked by the discovery of persistent emission from SGRs that is similar to the 
AXPs and bursting behaviour in the AXPs. See Woods & Thompson (2004) for a review. 

It is with the discovery of optical and infrared counterparts (Hulleman et al., 2000) 
that models for AXPs other than the magnetar became untenable. Hulleman et al. (2004) 
were the first to produce a spectral energy distribution for the AXP 4U 0142+61, showing 
an intriguing hint of a spectral feature: a sharp break in the optical. Furthermore, Kern & 
Martin (2002) discovered, optical pulsations with pulsed fraction of the order 25%, modulated 
at the X-ray period. 

Although in the magnetar model all the details are not yet ironed out, it does explain 
how such emission can arise from cyclotron emission by ions in the outer magnetosphere 
(Thompson, 2004, priv. comm.; Thompson, Lyutikov & Kulkarni, 2002). From the magnetar 
theory, it is not entirely clear whether one would expect the X-ray and infrared fluxes to be 
correlated. Such a relationship has been seen for the AXP IE 2259+586, which has recently 
been shown to have correlated X-ray and infrared fluxes following an X-ray bursting episode 
(Tarn et al., 2004). The return to the quiescent flux was found to occur on the same time- 
scales in both bands, and the flux ratio to remain roughly constant throughout the active 
episode. 

IE 1048.1-5937, is a 6.4s AXP in the field of the Carina Nebula, one of the two AXPs 
to have shown SGR-like outbursts to date and also the most noisy in terms of its timing char- 
acteristics (Kaspi et al., 2001; Gavriil, Kaspi & Woods, 2002). This, as well as its relatively 
hard spectrum, make it the most SGR-like of the AXPs (Woods & Thompson, 2004). A pos- 
sible infrared counterpart to IE 1048.1—5937 was observed by Wang & Chakrabarty (2002; 
henceforth WC02), who found the magnitudes: J = 21.7(3), H = 20.8(3) and K s = 19.4(3). 
Intriguingly, the source was not detected by Israel et al. (2002) down to limits of J ~ 23, 
H ~ 21.5 and K s ~ 20.7. Ergo the object is highly variable, and therefore of particular 
interest. 

Here we present further deep photometry in both the optical and the infrared, showing 
detections in three bands at levels far fainter than those of WC02, but consistent with Israel 
et al.'s limits. It may be that WC02 observed the object in a state similar to IE 2259+586 
in which its infrared flux was enhanced because of a preceding outburst. If so, this work 
probably shows the spectral energy distribution of the object in its quiescent phase. 

Below, we will first present imaging observations performed at the VLT and Magellan, 
and describe their analyses. Next, we describe the implications of these in three main areas: 
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the variability, spectral energy distribution, and X-ray to infrared flux ratio. We conclude 
with a brief discussion. 



2. Observations and Analysis 

Infrared imaging (J, H and K s bands) was performed of the field of AXP IE 1048.1—5937 
using ISAAC at the VLT. At Magellan, further imaging data were taken both in the optical 
(I, z' and a wide 'VR' filter) and infrared (J and K s ) using MagIC and PANIC, respectively. 
See Table 1 for a list of observations and the conditions for each, and Figure 1 for images of 
each of the detections. 



2.1. Astrometry 

The images below were referenced to the International Celestial Reference System 
through identifying stars in the USNO B1.0 catalogue on a short (30s) I-band image. Sixty- 
five stars were cross-identified, and a solution found with RMS deviations ~ 0.3" in each 
coordinate (after rejecting objects with residuals greater than 0.6", leaving 49 good matches), 
and thus an astrometric error of 0.3/\/49 = 0.04" in connecting to the USNO reference frame. 
This in turn has a systematic uncertainty of 0.2" (Monet et al. 2003) in connecting to the 
ICRS. Thus the total uncertainty in astrometry is ~ 0.2" in each coordinate. 

The uncertainties in connecting the coordinates of the short I-band exposure to the rest 
of the frames are negligible in comparison with the above errors. The 0.6" error radius, at 
90% confidence, in the Chandra position of IE 1048.1-5937 (Zombeck et al., 1995; WC02) 
combined in quadrature with the astrometric error above, corresponds to an error of ~ 0.74" 
on any of the images in Figure 1. 

2.2. VLT/ISAAC 

On the night of 23 April 2003, infrared images of the field of IE 1048.1—5937 were taken 
with ISAAC (Infrared Spectrometer And Array Camera; Moorwood, 1998) on Antu, unit 
telescope 1 of the Very Large Telescope. We used the Long Wavelength (LW) arm with the 
Aladdin lkxlk InSb infrared array, with pixel size 0.148". Seeing conditions were excellent. 
See Table 1 for a list of integration times and seeing for each image. The source was detected 
in the long J- and K s -, but not in the shorter H-band observations. 
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Right Ascension (J2000) Right Ascension (J2000) 

Fig. 1. — Images of the field of IE 1048.1—5937, in all bands with detections: J-band (top 
left), I-band (top right) and the K s -band (bottom) centred on the position of the AXP. The 
left-hand images were taken with the VLT and the right-hand ones with Magellan. The 0.74" 
radius error circle is shown on each image. Note that the horizontal band-like structure in 
the Magellan K s image is due to a irreproducible "pattern noise" that existed on the PANIC 
instrument at the time (see text). 
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Table I. Observation log of IE 1048.1-5937. 



Date Telescope, Mid-exposure Band Integration Seeing a Magnitude vF v 





instrument 


(UT) 




time (s) 


(arcsec) 




(erg s 1 cm 2 ) 


24 April 2003 


VLT 


01:40 b 


K s 


2772 


0.34 


21.3(3) 


2.6(9) x 10~ 15 




ISAAC 


01:45 b 


H 


1848 


0.47 


>21.3 


< 5.6 x 10" 15 






02:15 


J 


3255 


0.43 


23.4(4) 


1.6(8) x 10" 15 


6 June 2003 


Magellan 


03:50 


I 


7200 


0.33 


26.2(4) 


3.0(13) x 10~ 16 




MagIC 


05:35 


VR 


2700 


0.45 


>26.0 C 


< 6.8 x 10~ 16 






06:30 


z' 


2700 


0.44 


>24.2 d 


< 3.0 x 10" 15 


7 June 2003 


Magellan 


03:30 


K s 


2300 


0.38 


>20.9 


< 3.8 x 10" 15 




PANIC 


06:00 


H 


3780 


0.34 


>20.8 


< 9.0 x 10" 15 






05:00 


K s 


1375 


0.30 


21.5(4) 


2.2(10) x 10" 15 



Note. — Limits are at 95% confidence levels. 
a Full width at half maximum. 

b The VLT images were taken alternating between the filters, hence the close mid-exposure times. 

c This is on a Vega-like magnitude scale, with F V (VR = 0) = 3.33(2) x 10~ 20 erg s _1 cm~ 2 Hz -1 (see text). 

d This magnitude is not on the Vega, but on the AB-system (see text). 
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Standard reduction was carried out in order to subtract the dark current, flat field 
(using median averages derived from the science frames) and combine the frames using 
IRAF. Photometry was performed using DAOPHOT II (Stetson, 1987). 

In order to calibrate the frames, aperture photometry was performed on short-exposure 
images of three standard stars obtained on the same night. The aperture size was set large 
enough to include the majority of the flux from the star, as the sky noise was negligible in 
comparison; this established the magnitude zero point for each of the frames. The airmass 
of these exposures was in every case within 0.1 of the science exposures, so an atmospheric 
extinction correction was not necessary. Colour terms were negligible. In order to find the 
magnitude of the stars in the science field, the aperture correction was calculated using the 
difference between the magnitudes produced by the allstar task and those derived using 
the same large aperture used for the standards for several bright, relatively isolated stars in 
the frame (after the removal of fainter neighbours). See Table 2 for a list of the magnitudes 
of field stars around the AXP, with the numbering scheme as in WC02. 

The errors on measurements, and also the limiting magnitudes where the source was 
not detected, were determined empirically. A routine was run which inserted stars of known 
magnitude (scaled from the PSF) into areas of the image free of stars, and then the new 
image was passed through the same analysis procedure as for the science frame and the 
standard deviation of the measured magnitudes was calculated for the inserted stars. The 
latter is used as the measurement error for that known magnitude. It incorporates the fact 
that fainter stars will be spread over a wider magnitude range as measured, and that near 
the magnitude limit the faintest stars will not be measured at all. The method gets around 
any assumptions of the flat field accuracy or profile errors, which are assumed constant 
parameters in the error model used by DAOPHOT. Whilst such parameters could be fitted 
for each image in turn, and then used to determine the errors, we found that in every case 
the function of error versus magnitude found by DAOPHOT deviated both from our standard 
deviations and from photon noise whatever the choice of flat-field and profile error. This 
turned out to be significant at the very bright end and near the magnitude limit. 

We find, in the manner described above, the following magnitudes for the star XI, the 
AXP counterpart: H > 21.3, J = 23.4(4), K s = 21.3(3). The limit in H corresponds to 95% 
confidence. 



-7- 



2.3. Magellan/PANIC 

On the night of 6 June 2003, PANIC (Persson's Auxilliary Nasmyth Infrared Camera; 
Martini et al., 2004), the lkxlk infrared imaging array with 0.125" pixels on the Magellan 
Clay Telescope 1 , was used to acquire images of the field of IE 1048.1—5937 in the H- and 
K s -bands under good conditions. Two separate imaging runs were performed in the K s -band 
through the night, resulting in two images for analysis. The AXP was detected in only one 
of the K s frames and not at all in H; see table 1. 

There was a roughly periodic pattern noise in the PANIC detector at this time. This 
noise was not reproducible, and its effect can be seen in our final processed images in Figure 
1 as a band-like variation in the brightness of the sky. 

The images were processed and analysed as above, unsing three standard stars. We find 
the following magnitudes for the AXP: K s = 21.5(4), K s > 20.9, and H > 20.8. 

2.4. Magellan/MagIC 

On 5 June 2003, MagIC (The Raymond and Beverly Sackler Magellan Imaging Camera 2 ; 
Shectman & Johns, 2003), a 2kx2k CCD with quad readout amplifiers and 0.069" pixels 
on the Clay Telescope, Las Campanas was used to image IE 1048.1—5937 in the I- and 
z'-bands. A set of images were also taken using the custom 'VR' filter installed on MagIC 
(a wide, roughly rectangular pass band covering much of V and R; see Jewitt, Luu & Chen, 
1996). For the I band, about 90% of the integration time was at excellent seeing conditions of 
< 0.4", whereas the average seeing was 0.44" in VR and 0.46" in z'; see Table 1. The I-band 
images with the brightest sky and worst seeing were excluded from the stacking process, 
leaving a total of integration time of 7200s. 

Again, the frames were processed and combined using IRAF-DA0PH0T II, with bias sub- 
traction and trimming performed by the ccdmagic task provided by the observatory and flat 
fields again derived from the science frames (this proved more successful than screen flats). 
Magnitudes were calibrated using several standard stars and errors calibrated empirically as 
before (following the ISAAC experience above). 

Seven photometric standard stars for this run were taken from the E5 field (Stetson, 
2000), at an airmass close to the mean airmass of the science exposure. The I-band zero 



1 see http : //www . ociw . edu/lco/magellan/ instrument s/PANIC/ panic/ 
2 see http : //occult . mit . edu/ instrumentation/magic/ 
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point was a simple average over the values inferred from the well exposed stars on the frame. 

The source was not detected in the VR-band. For calibration of this non-standard filter, 
the E5 standard field was also imaged with the VR filter, and also the standard V and R 
filters to check for consistency. The VR zero point and effective wavelength were deduced in 
the following manner: the standard stars were assumed to have power-law spectra over their 
V. . . R range, so that the VR magnitude would be a simple interpolation, VR = (l—b)V+bR, 
where the parameter b is to be found. Next, VR was defined to be such that the colour-term 
coefficient is zero: VR — vr + zvr (where vr is the instrumental magnitude and zvr the 
zero point). In effect, this makes this particular VR filter the prototype for this magnitude 
scale. Finally, the value of b was found which minimised variance about an average value 
for the zero point. This minimum value was found to be b = 0.51(7), yielding the values: 
\ VR = 0.592(7)/zm, z VR = 27.9(4) (for le-/s), and F V (VR = 0) = 3.33(4) x 10~ 20 erg s" 1 
cm -2 Hz -1 . 

The z' band is on the magnitude AB system, and is based on a different set of standard 
stars for its basic calibration. In order to find the zero point, well-known empirical transfor- 
mations between z' and the VRI bands were used (Smith et al., 2001) to find the magnitudes 
of each of the standard stars. These were compared with the measured magnitudes of the 
standard field with the z' filter. 

The AXP counterpart was found to have magnitudes I = 26.2(4), VR > 26.0 and 
z' > 24.2. 



3. Results and Implications 

Figure 2 shows a colour-colour diagram for stars in the field of IE 1048.1—5937. The 
AXP is clearly offset from the bulk of the stars. Its J-K s = 2.1(5) colour is similar to that 
found by WC02 (who found J-K= 2.3(4)), despite the fact that it is ~ 2mag fainter, and 
to that of the brightest AXP, 4U 0142+61. Given the variability and peculiar colours, we 
believe there is no longer any doubt that the source, labelled XI by WC02, lying at the 
centre of the error circle derived from Chandra is indeed the infrared/optical counterpart to 
the AXP. 

We now discuss the inferred spectral energy distribution, X-ray to infrared flux ratio 
and variability. The spectral energy distribution and X-ray to infrared flux ratio of the AXP 
depend on our assumption for the amount of reddening to the source. Our knowledge on 
this matter is discussed in Section 3.2. 
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Table 2. Photometry for selected stars in the field of IE 0148.1-5937. 



Star id a 


I 


J 


H 


K s 


XI 


26.2(4) 


23.4(4) 


>21.5 


21.3(3) 


X2 


23.09(2) 


20.44(5) 


19.51(5) 


19.21(4) 


X3 


24.85(14) 


21.7(2) 


20.64(15) 


20.00(8) 


X4 


23.54(4) 


21.0(1) 


20.12(9) 


19.80(6) 


X5 


18.51(2) 


16.83(2) 


16.34(3) 


16.28(2) 


X6 


22.24(2) 


19.66(4) 


18.67(4) 


18.46(2) 


X7 


22.59(3) 


20.21(5) 


19.35(5) 


19.11(3) 


X8 


20.61(2) 


18.63(2) 


18.05(3) 


17.95(2) 


f 


22.95(3) 


20.46(5) 


19.61(5) 


19.24(3) 


A 


15.59(2) 


14.68(2) 


14.41(3) 


14.46(2) 


B 


16.47(2) 


14.97(2) 


14.65(3) 


14.59(2) 


C b 


17.42(2) 


15.96(2) 


15.52(3) 


15.54(2) 


D b 


17.11(2) 


14.57(2) 


13.69(3) 


13.46(2) 



Note. — I-band magnitudes are from the Magellan 
data, the rest from the VLT. 

a As labelled in WC02. Figure 1 shows the locations 
of some of these sources within the field. 

b There appears to have been an error in WC02, 
where the IR magnitudes for Stars C and D were 
switched (Wang, Z., 2005, priv. comm.) 
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Fig. 2.— Colour-colour plot of stars in the field of IE 1048.1-5937. The AXP (star XI) 
is labelled with an open square and the arrow shows the direction of increasing reddening. 
The outlier in the top-right is probably a genuine, highly reddened background star and 
a 3o" deviation; it lies far from the error circle. 
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Fig. 3. — Photometry for IE 0148.1—5937, as measured. The filled circles 
represent this work (where the H-, z'- and VR-band points are the respec- 
tive best limits, and for K s only the more accurate measurement, from the 
VLT, is shown); open triangles are those of WC02 and open squares the limits 
established by Israel et al. (2002). 
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3.1. Variability 

As has already been noted by Israel et al. (2002), IE 1048.1—5937 has shown large 
variability in its infrared flux. The Kg-band magnitudes presented here are much fainter 
than those given in WC02, but consistent with the limits found by Israel et al. This is 
illustrated in Figure 3. 

IE 1048.1—5937 is currently part of a regular monitoring programme using the Rossi 
X-ray Timing Explorer (RXTE), and has been observed many times over the last two years, 
at particularly closely spaced intervals during the time of the observations above (Gavriil 
& Kaspi, 2004). Marked changes are seen during this period in pulsed flux 3 and spin-down 
rate (z>). The hardness ratio shows no clear variations during this time. Figure 4 shows the 
pulsed flux and spin-down rate, with K-band photometry points as a comparison. One sees 
a hint of anti-correlation between the X-ray and infrared fluxes, and no obvious relationship 
between the infrared flux and spin-down rate. As Gavriil & Kaspi (2004) stated, X-ray flux 
and spin-down rate are clearly not related. 

For another source also in the RXTE monitoring programme, IE 2259+586, recent 
results indicate that X-ray and infrared fluxes were correlated in the period following multiple 
X-ray bursts and increased activity (Tam et al., 2004) with both spectral bands showing 
increased persistent fluxes, which decayed towards the pre-burst values on the same time- 
scale. This suggests that the mechanism at work in IE 2259+586 is different from that 
of IE 1048.1—5937, possibly because of the bursting activity (note that bursting activity 
cannot be excluded by the RXTE observations of IE 1048.1-5937). 

3.2. Distance and Reddening 

The intrinsic spectrum that we derive for the AXPs depends strongly on the value that 
we assume for the interstellar reddening along the line of sight. Since the reddening to all of 
the magnetars in the plane of the Galaxy have nominal reddening Ay > 5 (which includes 
all of the AXPs, apart from the recently discovered CXOU J010043. 1-721134, Majid et al., 
2004), the extinction correction is huge and has a profound impact on the spectrum derived. 
As an example, the inferred optical spectrum of 4U 0142+61 changes from Raleigh- Jeans 
like to flat (in vF u ) for a change Ay = 5 to Ay = 2 (Hulleman et al., 2004). It is therefore 
important to consider this as part of the analysis. 



3 Due to nearby strong X-ray sources and the non-imaging nature of the RXTE observations, only the 
pulsed flux could be measured 
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Fig. 4. — Spin-down rate, v (top panel, crosses), K s -band apparent magnitude (top panel, 
filled squares) and pulsed flux (lower panel) as a function of time. The spin-down and flux 
data are taken from Gavriil & Kaspi (2004). The wide point centred at (MJD = 51650, z> = 
—9.15 x 10~ 13 ) represents a period in 2000 throughout which phase-connected timing was 
possible (Kaspi et al., 2001). At MJD~52400, the observing strategy was changed to enable 
spin-down measurements without long time-baseline phase coherence. 
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Usually, the model for the intrinsic X-ray spectrum is used to infer the amount of 
extincting material (given as the Hydrogen column, Njj), and assuming the relationship 
derived by Predehl & Schmitt (1995) to calculate the reddening due to dust. 

For IE 1048.1—5937, Tiengo et al. (2002) found an apparently featureless X-ray spec- 
trum well-fitted by an absorbed power law [yF v ~ v~ a , a ~ 2 — 3) plus blackbody (kT ~ 0.64 
keV). This model for the X-ray spectrum yields an extremely high value for the extincting 
column to the source of Nh = 1.04(8) x 10 22 cm -2 , comparable to values of total extinction 
seen through many lines of sight through the galactic plane (but see also below). In this 
way, values typically of order Ay ~ 5.8 have been arrived at (WC02) for IE 1048.1—5937. 
Comparing this estimate to the typical reddening to stars in the Carina nebula (Ay ~ 2), 
Seward et al. (1986) argued that the AXP must lie behind the Carina Nebula, and thus 
placed a lower limit on its distance of d > 2.8 kpc. 

There exist other independent methods for estimating the run of reddening with distance 
in any given line of sight. For our field, we found the following relevant measurements. 

Neckel, Klare & Sarcander (1980) measured the average function of reddening with 
distance, Ay(d), in many galactic-plane fields including one containing this source. They 
find that as a function of distance, the visual extinction is constant between 1.5kpc - 6kpc at 
Ay ~ 1. It should be stressed, however, that this value is an average over their field (size of 
order ~square degrees), and one could argue that measurement biases would favour a lower 
value of extinction at any given distance in an area of highly variable dust content. 

Carrero et al. (2004) performed colour-magnitude studies of three nearby open clusters 
in the Carina complex (approximately 1° from the region of interest). For the clusters Trum- 
pler 16 and Trumpler 14, they find that Ay =2.0(1) and 1.8(7) and distances of d =3.9(5) 
and 2.5(3) respectively. This is in a region where the dust density is visibly greater than 
the surrounding sky; for IE 1048.1—5937 to have a greater extinction, this suggests a more 
constraining lower bound on its distance. It should be noted however, that the authors of 
this study stress the variability of reddening from star to star in the clusters. 

Finally, Merehgetti, Caravea & Bignami (1992) find that one of the brighter stars in the 
field is probably very highly redenned (Ay > 7), based on two interstellar absorption bands 
(4428 and 5778A) in its spectrum. The distance to this star is unknown. 

To summarise, there is reason to question our assumption that the value of N H gained 
through the fitting of the X-ray spectrum gives an appropriate value of reddening; such a 
suggestion has also arisen in the case of the brightest of the AXPs, 4U 0142+61. Hulleman, 
van Kerkwijk & Kulkarni (2004) suggest that values in the range Ay = 2 ... 6 are possible. 
A similar uncertainty was also suggested in the case of AXP IE 2259+586 (from the mea- 
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surements of optical filaments in its associated supernova remnant, Fesen & Hurford, 1995). 
As we see below, the value we use to de-redden flux measurements will alter the shape of 
the resultant spectrum drastically. 

3.3. Spectral Energy Distribution 

Figure 5 shows the inferred spectral energy distribution (SED), de-redenned using Ay = 
5.8 (the value derived from the X-ray spectral fit). As can be seen, a line of vF v ~ const 
would be consistent with these points (the H-and VR-band limits are consistent with such a 
hypothesis). 

Also on Figure 5 is plotted the SED of 4U 0142+61 (Hulleman, van Kerkwijk & Kulkarni, 
2004), showing the optical and infrared regions for comparison, de-redenned with the nominal 
value for that source of Ay = 5.1. The inferred SEDs of the two AXPs are clearly inconsistent 
(see particularly the I-band point). 

It is worth noting here, following from Section 3.2, that the reddening to both these 
objects is rather uncertain. With Ay = 5.8, IE 1048.1-5937 has as SED with vF v = u a , 
a = 0.2(4). This changes to a = —1.3(4) for Ay ~ 2. Hulleman, van Kerkwijk & Kulkarni 
(2004) suggest that for 4U 0142+61, values of reddening as low as Ay ~ 2 are possible, 
based on the run of reddening with distance along the line of sight. Were one to take 
Ay ~ 2 for 4U 0142+61, its SED would look much flatter and not unlike that which is 
shown for IE 1048.1—5937 with Ay = 5.8. Thus it could indeed be that the two have similar 
intrinsic spectra but that the value assumed for reddening is wrong for either or both of the 
sources. 

If we were to assume that the intrinsic spectra of the two AXPs IE 1048.1—5937 and 
4U 0142+61 are the same, then this allows us to calculate the relative reddening between 
the two objects, AA V . This value enables a comparison between the two objects without 
knowledge of the specific reddening to either one. By comparing the I — K colours of the 
two objects we find AE(I-K) = (I - K) 1048 - (I - K) Q142 =1.2(5), implying AAy = 2.5(5). 
Thus, if the intrinsic spectra are indeed the same, IE 1048.1—5937 would have V = 28.2(5) 
and R = 26.9. These predictions can be tested. 

3.4. X-ray to Infrared Flux Ratio 

In Table 3 we list the currently known X-ray to infrared flux ratios of the AXPs, corrected 
for extinction using the nominal values of Nh derived from fits to the X-ray spectra. Note 
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Fig. 5. — Spectral energy distribution for IE 0148.1—5937 (circles), both 
as measured (open) and de-reddened with Ay = 5.8 (filled). Also shown is 
the spectral energy distribution for the brighter AXP 4U 0142+61 (squares), 
multiplied by a factor of 100 for clarity and de-redenned with A v = 5.1 
(filled squares). Only detections are shown; notice that the B-band point 
for 4U 0142+61 is a marginal detection. The data for 4U 0142+61 is taken 
from Hulleman, van Kerkwijk & Kulkarni (2004) and Israel et al.(2004). 
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that neither the 2. . . lOkeV range nor the K-band suffer considerable extinction, so these 
numbers are not heavily dependent on the values of reddening chosen. However, some of the 
sources have shown large variability in both their X-ray and infrared fluxes, and the numbers 
here represent values only at some particular time. 

Given the variability, the flux ratios are remarkably consistent, except for that of 
1RXS J170849— 400910. Interestingly, the counterpart for the latter is the least secure 
amongst the AXPs, as its association is based on astrometry and possible peculiar colours 
only (Israel et al., 2003). In fact, its inconsistency with the other sources in this table sug- 
gests that perhaps the other, fainter source within the Chandra error circle (named Star 'B' 
by these authors) could be the counterpart. The relative contributions of the sources to the 
J-band flux is unknown, and might account for the apparent strange colours of Star A. 

IE 1048.1—5937 has varied in both its X-ray and infrared emission (see Section 3.1) 
by over a factor of 2, so the range which the flux ratio could take, when measured at some 
particular epoch is large. From Figure 4, we see that the ratio of pulsed X-ray flux to infrared 
flux has varied from 8 x 10 12 . . . 1 x 10 14 (in units of ctsPCU^ 1 / ergcur 2 ) . On the one hand, 
then, the X-ray to infrared flux ratios of AXPs are similar, and for IE 2259+586 the X-ray 
and infrared fluxes were correlated following an outburst. On the other hand, the infrared 
and X-ray fluxes of IE 1048.1—5937 seem, if anything, to be anti-correlated. These two 
statements appear hard to reconcile with one-another. 

4. Discussion 

In this work, we have detected the likely quiescent counterpart to the AXP IE 1048.1—5937. 
To summarise our findings, we have found large (~2mag) variations in the infrared bright- 
ness of IE 1048.1—5937 which suggest a possible anti-correlation with X-ray flux, but no 
clear relationship with spin-down rate. This variability is the largest that has been seen in 
the AXPs, but variability is not unique to this source. The anti-correlation between infrared 
and X-ray flux is opposite to the behaviour of IE 2259+586, where a positive correlation 
was found (Tarn et al., 2004). 

The infrared spectrum has remained of a consistent shape with WC02. We find that the 
inferred spectral energy distribution of IE 1048.1—5937 is inconsistent with that inferred for 
4U 0142+61, but that they can be made consistent with an appropriate choice of relative 
reddening (for which there may be justification). We find that although both the infrared and 
X-ray for this source vary a large amount and are not correlated, that the X-ray to infrared 
flux ratio is consistent with the other AXPs. Perhaps this challenges our understanding 
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Table 3. Comparison of the un-absorbed X-ray flux and infrared flux for all Anomalous 

X-ray Pulsars with infrared counterparts. 



AXP 


X-ray flux a 






K-band flux d 


Fx /Fx 




(erg s^ 1 cm~ 2 ) 




(10 22 cm" 2 ) 


(erg s _1 cm" 2 ) 


/1000 


IE 1048.1-5937 


1.4 x 10- 11 


21.3 


1.0 


5.1 x 10~ 15 


2.7 


4U 0142+61 


8.3 x 10- 11 


20.1 


0.91 


1.4 x 10~ 14 


5.9 


1RXS Jl 70849-4009 10 


6.4 x 10~ n 


18.3 


1.4 


9.7 x 10~ 14 


0.66 


IE 2259+586 


2.0 x 10~ n 


21.7 


1.1 


3.5 x 10~ 15 


5.7 


XTE J1810— 197 


2.2 x 10~ n 


20.8 


1.1 


8.1 x 10~ 15 


2.7 



Note. - - Data from Woods & Thompson (2004) and references therein. The values are 
not thought to be within active periods. 

a Un- absorbed flux in the 2-10keV range. 

b Apparent K s magnitude. 

c Inferred from fits to the X-ray spectrum. 

d vF I/ , de-redenned assuming A v = N H /1.79 x 10 21 (Predehl & Schmitt, 1995). 
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of whether the AXPs do indeed have quiescent periods, and of the time-scales on which 
the magnetosphere reacts to internal magnetic changes. We now discuss what relationships 
might be expected from the magnetar model. 

The model proposed for how infrared and optical emission might arise from a magnetar 
is as follows. The decay of the internal magnetic field of the magnetar induces a twist in 
the magnetosphere of the magnetar. This large-scale twist is supported by global, persistent 
currents, which exist out to large radii. These currents power the non-thermal X-ray and 
gamma-ray luminosity, and may also produce infrared and optical emission, through ion 
cyclotron radiation at large radii (Thompson, Lyutikov & Kulkarni, 2002). 

In the case that the infrared radiation is indeed cyclotron emission by ions gyrating 
in the outer magnetosphere, the infrared flux and torque should be correlated because of 
the following argument. First, the number of ions travelling to large enough radii to have 
Landau transition energies in the optical/infrared range is a fixed fraction of the total electric 
current. The total current in the magnetosphere defines the amount of twist in the global 
magnetic field structure (that is, the departure from a pure dipolar field), which in turn 
affects the rate at which the field strength falls off with radius. Since the torque is caused 
by the field strength at the light cylinder, it must then be correlated to the ion population 
at large radii and so to the infrared flux (Thompson, Lyutikov & Kulkarni, 2002). 

Whether the infrared and X-ray flux are expected to be correlated depends on the 
details of their respective emission mechanisms. Changes in X-ray flux would be possible, 
for example, from internal heating, which would not directly affect the infrared emission. 

In order to constrain the emission mechanism for the infrared flux, it will be important 
to establish whether the anti-correlation between infrared and X-ray flux is indeed real, and 
whether the IR flux to torque correlation can really be ruled out. It would also be interesting 
to look for a correlation with the flux in the strong, hard X-ray tail of the spectrum that was 
recently discovered by Kuiper at al. (2004). Furthermore, it would be worthwhile to seek 
further spectral signatures, and to find or at least constrain the reddening and distance to 
the various sources. 
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